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ABSTRACT 

Using a propagation calculation with energy dependent 
parameters, including the depletion of short pathlengths, 
and incorporating experimental nuclear excitation functions, 
the variation of the °Ber Be ratio with the matter 
densities in two nested confinement regions is investigated. 
It is shown that there is no unique correspondence Between a 
^Be/%e measurement at low energy and the density of matter 
in the galaxy. ^°Be/°Be measurements at both low and high 
energy are needed to fully specify the matter densities. 


1. Introduction In recent years considerable progress has been made 
in determining the energy dependence of the pathlength distribution 
(PLD) for cosmic ray propagation at low energies in the interstellar 
medium (Garcia-Munoz et al., 1985)- However, previous analyses of 
radioactive cosmic ray isotopes have been done in the context of the 
simple "leaky box" model, which does not explain, simultaneously, the 
energy dependence of the measurements of both light and heavy secondary 
to primary elements (such as B/C and V/Fe or Sc/^g). In this report, we 
present results of propagation calculations for Be/^Be using 
experimental nuclear excitation functions and energy dependent 
parameters in the PLD which give results in agreement with both B/C and 
sub-Fe/Fe ratios from 100 MeV/n to 30 GeV/n. 

2. The Propagation Calculations The propagation code employs the 

weighted slab technique to calculate the abundances of 96 stable, long- 

lived or electron capture isotopes from **He to Radioactive 

decay (B , B _ , electron capture) is treated explicitly; the effects of 

ionization energy loss are included; and energy dependent total 
inelastic cross sections, based upon compiled experimental data, are 

employed. The partial fragmentation cross sections are based upon semi- 

empirical formulae of Silberberg and Tsao (1973), modified using 

available experimental data. Cosmic ray source abundances are from 

Garcia-Munoz and Simpson (1979), and the ratio of H to He in the 
interstellar medium is from Cameron (1981). The majority of the cosmic 
ray data used here were collected during the last period of minimum 
solar modulation (1974-1979), and the calculations include modulation 
with an adiabatic deceleration parameter <(> = 490 MV (Evenson et al., 
1984). 
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The propagation 
code was run in two 
modes: single and double 
runs. In a single run, 
the abundances of all 
species were calculated 
for a series of slabs 
and then weighted by a 
single, energy dependent 
pathlength probability 
distribution. Such a 
PLD is illustrated in 
Fig, 1 where the inset 
(B) shows the overall 
shape of the Double Exponential 



Kinetic Energy (MeV/nucleon) 

Fig. 1: The energy dependent DE PLD. 

(DE) form, composed of two exponentials 


with means X-, and X 2 , where X 2 represents the depeletion of short 
pathlengths tGarcia-Munoz et al., 1984). The energy dependences of X-j^ 
and X 2 , shown on part A along with the mean pathlength <X>, are 
determined by fitting the B/C and sub-Fe/Fe data over the energy 
range 0.1 - 30 GeV/nucleon (Guzik and Wefel, 1984). The DE PLD is 
qualitatively similar to earlier models (Garcia-Munoz and Simpson, 
1970), e.g. the PLD in the "nested leaky-box" model (Cowsik and Wilson, 
1975), but is quantitatively different since our DE PLD includes energy 


dependent parameters. 


Fig, 2 shows a comparison of the calculated Be/%e ratio for 
different PLD’s, for a constant density of 0,5 atoms /cm. The dashed 
line is for a simple energy independent exponential which is 
characteristic of the "leaky box" model. Note that when energy 

the PLD is included (lower solid line) the 
Be/^Be is modified. However, neither of these 
PLD* s fit both the B/C and sub-Fe/Fe ratios. Such a fit requires an 
energy dependent depletion of short pathlengths (a truncated PLD, see 
Garcia-Munoz et al., 1984) and results for two types of truncation 


dependence of the mean of 
shape of the predicted 


(yielding essentially identical 
results) are shown on Fig. 2 
(upper solid and dot-dashed 
curves). Note that above 
several GeV/nucleon the curves 
converge, as expected, since at 
higher energies the truncation 
becomes negligable (c.f. X 2 on 
Fig. l). 

In analogy with the "nested 
leaky box" model, the two 
components of the energy dependent 
DE PLD can be associated with 
confinement in an inner region of 
density p (i.e. around the 
sources) nested within an outer 
region of density P 0T1TT1 (i.e. the 
Galaxy) (Guzik and Wefel, 1984). 



KINETIC ENERGY (MeV/nucleon) 

Fig. 2: “^Be/^Be results for 

different PLD's. 
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Fig. 3: ^Be/^Be ratio as a 

function of density. 


Double runs were used to study this 
configuration of confinement regions. 

In the first step, the propagation code 
is run with an exponential PLD whose 
mean follows the curve given by ^ on 
Fig. 1. For this step the matter 
density has the value p . Next, the 
results of step 1 are used as the 
source for a second calcualtion (with 
matter density p Q ) whose PLD mean 
follows the curve on Fig. 1. This 
two step method allows us to treat 
separately the densities in the inner 
and outer regions, which, physically, 
could be quite different. 



KINETIC ENERGY (MeV/nucleon) 


3. Results 


Figure 3 shows 


for different matter densities 
single runs using the DE PLD with 
P™ = P, 


results 

from 


Fig. k: Results of the two- 

step analysis compared to 
experimental data. 


"IN ^OUT* curves converge at 

high energy due to relativistic effects, 

and for low energy, the survival of u Be is directly dependent upon the 
matter density. The low energy satellite experiments give a mean value 
°Be/^Be = 0.13 ± 0.03 at ~100 MeV/n which implies from Fig. 3, a 
matter density of 0.23 ± 0.06 atoms/cm , consistent with previous 
results which employed single density analysis. 


Fig. k shows the result of the two-step calculations compared to 
experimental data (o, • - Garcia-Munoz et al., 1977, 1981; 

□ - Wiedenbeck and Greiner, 1980; V - Webber et al., 1977; ▼ -Webber and 
Kish, 1979; A - Hagen et al., 1977)* The top panel shows the 
calculated 10 Be/^Be ratio for fixed with p^ allowed to vary from 

0.05 to lCr atom/cm . Different values for p^ spread the predicted 
U BerBe ratio at low energy by about a factor of two. The curves 
converge above 1 GeV/n since the effect of the inner region decreases 
with increasing energy. In this model a measurement of the u Be/^Be 
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ratio does not uniquely determine the density in the confinement region. 

The measured 10 Be/^Be ratio at low energy can provide bounds on the 
density in the outer region. This is illustrated in the lower panel of 
Fig. 4 where curves A and B show the total spread in the predicted ratio 
for P = 0.2 atoms /cm. Assuming a large value for P™ allows a 
jailer value of P^^ to be used to reproduce a given value of the 
UBe/'UBe ratio, as shown by curve C which forms a lower bound to the 
experimental data. Conversely, assuming a minimum value for P^, a 
value of P 0OT “ 0.5 atom/cnr give an upper limit to the low energy data. 
Thus, the available measurements constrain P to the range 0.05 - 0.50 
atoms /cm , below the average density of ~1 atom/cnr for the galactic 
disk and above the density of 0.01 atoms /cnr of a galactic halo. 

4. Conclusions The cosmic ray PLD, including an energy dependent 
depletion of short pathlengths, can be represented as two "nested” 
^gnfinement regions having different matter densities. The available 
Ber Be measurements do not uniquely determine the density in the outer 
region, but jj-j^ey do limit the allowed values to a range of 0.05 and 0.50 
atoms/cm . u Be/ y Be data at high energy (currently unavailable) 
combined with existing measurements at low energy can determine the 
matter density in both volumes. 
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